Approaching Mass Cytometry Translational Studies by Experimental and Data Curation Settings.

Clinical studies are conducted to better understand the pathological mechanism of diseases and to find biomarkers associated with disease activity, drug response, or outcome prediction. Mass cytometry (MC) is a high-throughput single-cell technology that measures hundreds of cells per second with more than 40 markers per cell. Thus, it is a suitable tool for immune monitoring and biomarker discovery studies. Working in translational and clinical settings requires a careful experimental design to minimize, monitor, and correct the variations introduced during sample collection, preparation, acquisition, and analysis. In this review, we will focus on these important aspects of MC-related experiments and data curation in the context of translational clinical research projects.

A genome-wide association study for survival from a multi-centre European study identified variants associated with COVID-19 risk of death.

The clinical manifestations of SARS-CoV-2 infection vary widely among patients, from asymptomatic to life-threatening. Host genetics is one of the factors that contributes to this variability as previously reported by the COVID-19 Host Genetics Initiative (HGI), which identified sixteen loci associated with COVID-19 severity. Herein, we investigated the genetic determinants of COVID-19 mortality, by performing a case-only genome-wide survival analysis, 60 days after infection, of 3904 COVID-19 patients from the GEN-COVID and other European series (EGAS00001005304 study of the COVID-19 HGI). Using imputed genotype data, we carried out a survival analysis using the Cox model adjusted for age, age2, sex, series, time of infection, and the first ten principal components. We observed a genome-wide significant (P-value < 5.0 × 10-8) association of the rs117011822 variant, on chromosome 11, of rs7208524 on chromosome 17, approaching the genome-wide threshold (P-value = 5.19 × 10-8). A total of 113 variants were associated with survival at P-value < 1.0 × 10-5 and most of them regulated the expression of genes involved in immune response (e.g., CD300 and KLR genes), or in lung repair and function (e.g., FGF19 and CDH13). Overall, our results suggest that germline variants may modulate COVID-19 risk of death, possibly through the regulation of gene expression in immune response and lung function pathways.

Molecular characterisation of lupus low disease activity state (LLDAS) and DORIS remission by whole-blood transcriptome-based pathways in a pan-European systemic lupus erythematosus cohort.

OBJECTIVES: To unveil biological milieus underlying low disease activity (LDA) and remission versus active systemic lupus erythematosus (SLE). METHODS: We determined differentially expressed pathways (DEPs) in SLE patients from the PRECISESADS project (NTC02890121) stratified into patients fulfilling and not fulfilling the criteria of (1) Lupus LDA State (LLDAS), (2) Definitions of Remission in SLE remission, and (3) LLDAS exclusive of remission. RESULTS: We analysed data from 321 patients; 40.8% were in LLDAS, and 17.4% in DORIS remission. After exclusion of patients in remission, 28.3% were in LLDAS. Overall, 604 pathways differed significantly in LLDAS versus non-LLDAS patients with an false-discovery rate-corrected p (q)<0.05 and a robust effect size (dr)≥0.36. Accordingly, 288 pathways differed significantly between DORIS remitters and non-remitters (q<0.05 and dr≥0.36). DEPs yielded distinct molecular clusters characterised by differential serological, musculoskeletal, and renal activity. Analysis of partially overlapping samples showed no DEPs between LLDAS and DORIS remission. Drug repurposing potentiality for treating SLE was unveiled, as were important pathways underlying active SLE whose modulation could aid attainment of LLDAS/remission, including toll-like receptor (TLR) cascades, Bruton tyrosine kinase (BTK) activity, the cytotoxic T lymphocyte antigen 4 (CTLA-4)-related inhibitory signalling, and the nucleotide-binding oligomerization domain leucine-rich repeat-containing protein 3 (NLRP3) inflammasome pathway. CONCLUSIONS: We demonstrated for the first time molecular signalling pathways distinguishing LLDAS/remission from active SLE. LLDAS/remission was associated with reversal of biological processes related to SLE pathogenesis and specific clinical manifestations. DEP clustering by remission better grouped patients compared with LLDAS, substantiating remission as the ultimate treatment goal in SLE; however, the lack of substantial pathway differentiation between the two states justifies LLDAS as an acceptable goal from a biological perspective.

Type 1 interferons: A target for immune-mediated inflammatory diseases (IMIDs).

The improved understanding of the molecular basis of innate immunity have led to the identification of type I interferons (IFNs), particularly IFN-α, as central mediators in the pathogenesis of several Immune-mediated inflammatory diseases (IMIDs) such as systemic lupus erythematosus (SLE), systemic sclerosis, inflammatory myositis and Sjögren's syndrome. Here, we review the main data regarding the opportunity to target type I IFNs for the treatment of IMIDs. Type I IFNs and their downstream pathways can be targeted pharmacologically in several manners. One approach is to use monoclonal antibodies against IFNs or the IFN-receptors (IFNARs, such as with anifrolumab). The downstream signaling pathways of type I IFNs also contain several targets of interest in IMIDs, such as JAK1 and Tyk2. Of these, anifrolumab is licensed and JAK1/Tyk2 inhibitors are in phase III trials in SLE. Targeting IFN-Is for the treatment of SLE is already a reality and in the near future may prove useful in other IMIDs. IFN assays will find a role in routine clinical practice for the care of IMIDs as further validation work is completed and a greater range of targeted therapies becomes available.

Association of Combined Anti-Ro52/TRIM21 and Anti-Ro60/SSA Antibodies With Increased Sjögren Disease Severity Through Interferon Pathway Activation.

OBJECTIVE: The biologic diagnosis of primary Sjögren disease (SjD) mainly relies on anti-Ro60/SSA antibodies, whereas the significance of anti-Ro52/TRIM21 antibodies currently remains unclear. The aim of this study was to characterize the clinical, serological, biologic, transcriptomic, and interferon profiles of patients with SjD according to their anti-Ro52/TRIM21 antibody status. METHODS: Patients with SjD from the European PRECISESADS (n = 376) and the Brittany Diagnostic Suspicion of primitive Sjögren's Syndrome (DIApSS); (n = 146) cohorts were divided into four groups: double negative (Ro52-/Ro60-), isolated anti-Ro52/TRIM21 positive (Ro52+), isolated anti-Ro60/SSA positive (Ro60+), and double-positive (Ro52+/Ro60+) patients. Clinical information; EULAR Sjögren Syndrome Disease Activity Index, a score representing systemic activity; and biologic markers associated with disease severity were evaluated. Transcriptome data obtained from whole blood by RNA sequencing and type I and II interferon signatures were analyzed for PRECISESADS patients. RESULTS: In the DIApSS cohort, Ro52+/Ro60+ patients showed significantly more parotidomegaly (33.3% vs 0%-11%) along with higher ß2-microglobulin (P = 0.0002), total immunoglobulin (P < 0.0001), and erythrocyte sedimentation rate levels (P = 0.002) as well as rheumatoid factor (RF) positivity (66.2% vs 20.8%-25%) compared to other groups. The PRECISESADS cohort corroborated these observations, with increased arthritis (P = 0.046), inflammation (P = 0.005), hypergammaglobulinemia (P < 0.0001), positive RF (P < 0.0001), leukopenia (P = 0.004), and lymphopenia (P = 0.009) in Ro52+/Ro60+ patients. Cumulative EULAR Sjögren Syndrome Disease Activity Index results further confirmed these disparities (P = 0.002). Transcriptome analysis linked anti-Ro52/TRIM21 antibody positivity to interferon pathway activation as an underlying cause for these clinical correlations. CONCLUSION: These results suggest that the combination of anti-Ro52/TRIM21 and anti-Ro60/SSA antibodies is associated with a clinical, biologic, and transcriptional profile linked to greater disease severity in SjD through the potentiation of the interferon pathway activation by anti-Ro52/TRIM21 antibodies.

Local Ancestry at the Major Histocompatibility Complex Region is Not a Major Contributor to Disease Heterogeneity in a Multiethnic Lupus Cohort.

OBJECTIVE: Systemic lupus erythematosus (SLE) is an autoimmune disease resulting in debilitating clinical manifestations that vary in severity by race and ethnicity with a disproportionate burden in African American, Mestizo, and Asian populations compared with populations of European descent. Differences in global and local genetic ancestry may shed light on the underlying mechanisms contributing to these disparities, including increased prevalence of lupus nephritis, younger age of symptom onset, and presence of autoantibodies. METHODS: A total of 1,139 European, African American, and Mestizos patients with SLE were genotyped using the Affymetrix LAT1 World array. Global ancestry proportions were estimated using ADMIXTURE, and local ancestry was estimated using RFMIXv2.0. We investigated associations between lupus nephritis, age at onset, and autoantibody status with both global and local ancestry proportions within the Major Histocompatibility Complex region. RESULTS: Our results showed small effect sizes that did not meet the threshold for statistical significance for global or local ancestry proportions in either African American or Mestizo patients with SLE who presented with the clinical manifestations of interest compared with those who did not. CONCLUSION: These findings suggest that local genetic ancestry within the Major Histocompatibility Complex region is not a major contributor to these SLE manifestations among patients with SLE from admixed populations.

Serum profiling identifies CCL8, CXCL13, and IL-1RA as markers of active disease in patients with systemic lupus erythematosus.

Introduction: Systemic lupus erythematosus (SLE) is a clinically heterogeneous disease that presents a challenge for clinicians. To identify potential biomarkers for diagnosis and disease activity in SLE, we investigated a selected yet broad panel of cytokines and autoantibodies in patients with SLE, healthy controls (HC), and patients with other autoimmune diseases (AIDs). Methods: Serum samples from 422 SLE patients, 546 HC, and 1223 other AIDs were analysed within the frame of the European PRECISESADS project (NTC02890121). Cytokine levels were determined using Luminex panels, and autoantibodies using different immunoassays. Results: Of the 83 cytokines analysed, 29 differed significantly between patients with SLE and HC. Specifically, CCL8, CXCL13, and IL-1RA levels were elevated in patients with active, but not inactive, SLE versus HC, as well as in patients with SLE versus other AIDs. The levels of these cytokines also correlated with SLE Disease Activity Index 2000 (SLEDAI-2K) scores, among five other cytokines. Overall, the occurrence of autoantibodies was similar across SLEDAI-2K organ domains, and the correlations between autoantibodies and activity in different organ domains were weak. Discussion: Our findings suggest that, upon validation, CCL8, CXCL13, and IL-1RA could serve as promising serum biomarkers of activity in SLE.

Gut epithelial barrier dysfunction in lupus triggers a differential humoral response against gut commensals.

Introduction: Systemic lupus erythematosus is an autoimmune disease with multisystemic involvement including intestinal inflammation. Lupus-associated intestinal inflammation may alter the mucosal barrier where millions of commensals have a dynamic and selective interaction with the host immune system. Here, we investigated the consequences of the intestinal inflammation in a TLR7-mediated lupus model. Methods: IgA humoral and cellular response in the gut was measured. The barrier function of the gut epithelial layer was characterised. Also, microbiota composition in the fecal matter was analysed as well as the systemic humoral response to differential commensals. Results: The lupus-associated intestinal inflammation modifies the IgA+ B cell response in the gut-associated lymphoid tissue in association with dysbiosis. Intestinal inflammation alters the tight junction protein distribution in the epithelial barrier, which correlated with increased permeability of the intestinal barrier and changes in the microbiota composition. This permeability resulted in a differential humoral response against intestinal commensals. Discussion: Lupus development can cause alterations in microbiota composition, allowing specific species to colonize only the lupus gut. Eventually, these alterations and the changes in gut permeability induced by intestinal inflammation could lead to bacterial translocation.

Lack of strong innate immune reactivity renders macrophages alone unable to control productive Varicella-Zoster Virus infection in an isogenic human iPSC-derived neuronal co-culture model.

With Varicella-Zoster Virus (VZV) being an exclusive human pathogen, human induced pluripotent stem cell (hiPSC)-derived neural cell culture models are an emerging tool to investigate VZV neuro-immune interactions. Using a compartmentalized hiPSC-derived neuronal model allowing axonal VZV infection, we previously demonstrated that paracrine interferon (IFN)-α2 signalling is required to activate a broad spectrum of interferon-stimulated genes able to counteract a productive VZV infection in hiPSC-neurons. In this new study, we now investigated whether innate immune signalling by VZV-challenged macrophages was able to orchestrate an antiviral immune response in VZV-infected hiPSC-neurons. In order to establish an isogenic hiPSC-neuron/hiPSC-macrophage co-culture model, hiPSC-macrophages were generated and characterised for phenotype, gene expression, cytokine production and phagocytic capacity. Even though immunological competence of hiPSC-macrophages was shown following stimulation with the poly(dA:dT) or treatment with IFN-α2, hiPSC-macrophages in co-culture with VZV-infected hiPSC-neurons were unable to mount an antiviral immune response capable of suppressing a productive neuronal VZV infection. Subsequently, a comprehensive RNA-Seq analysis confirmed the lack of strong immune responsiveness by hiPSC-neurons and hiPSC-macrophages upon, respectively, VZV infection or challenge. This may suggest the need of other cell types, like T-cells or other innate immune cells, to (co-)orchestrate an efficient antiviral immune response against VZV-infected neurons.

Distinct gene dysregulation patterns herald precision medicine potentiality in systemic lupus erythematosus.

OBJECTIVES: We aimed at investigating the whole-blood transcriptome, expression quantitative trait loci (eQTLs), and levels of selected serological markers in patients with SLE versus healthy controls (HC) to gain insight into pathogenesis and identify drug targets. METHODS: We analyzed differentially expressed genes (DEGs) and dysregulated gene modules in a cohort of 350 SLE patients and 497 HC from the European PRECISESADS project (NTC02890121), split into a discovery (60%) and a replication (40%) set. Replicated DEGs qualified for eQTL, pathway enrichment, regulatory network, and druggability analysis. For validation purposes, a separate gene module analysis was performed in an independent cohort (GSE88887). RESULTS: Analysis of 521 replicated DEGs identified multiple enriched interferon signaling pathways through Reactome. Gene module analysis yielded 18 replicated gene modules in SLE patients, including 11 gene modules that were validated in GSE88887. Three distinct gene module clusters were defined i.e., "interferon/plasma cells", "inflammation", and "lymphocyte signaling". Predominant downregulation of the lymphocyte signaling cluster denoted renal activity. By contrast, upregulation of interferon-related genes indicated hematological activity and vasculitis. Druggability analysis revealed several potential drugs interfering with dysregulated genes within the "interferon" and "PLK1 signaling events" modules. STAT1 was identified as the chief regulator in the most enriched signaling molecule network. Drugs annotated to 15 DEGs associated with cis-eQTLs included bortezomib for its ability to modulate CTSL activity. Belimumab was annotated to TNFSF13B (BAFF) and daratumumab was annotated to CD38 among the remaining replicated DEGs. CONCLUSIONS: Modulation of interferon, STAT1, PLK1, B and plasma cell signatures showed promise as viable approaches to treat SLE, pointing to their importance in SLE pathogenesis.

A Summary on the Genetics of Systemic Lupus Erythematosus, Rheumatoid Arthritis, Systemic Sclerosis, and Sjögren's Syndrome.

Systemic lupus erythematosus, systemic sclerosis, rheumatoid arthritis, and Sjögren's syndrome are four major autoimmune rheumatic diseases characterized by the presence of autoantibodies, caused by a dysregulation of the immune system that leads to a wide variety of clinical manifestations. These conditions present complex etiologies strongly influenced by multiple environmental and genetic factors. The human leukocyte antigen (HLA) region was the first locus identified to be associated and still represents the strongest susceptibility factor for each of these conditions, particularly the HLA class II genes, including DQA1, DQB1, and DRB1, but class I genes have also been associated. Over the last two decades, the genetic component of these disorders has been extensively investigated and hundreds of non-HLA risk genetic variants have been uncovered. Furthermore, it is widely accepted that autoimmune rheumatic diseases share molecular disease pathways, such as the interferon (IFN) type I pathways, which are reflected in a common genetic background. Some examples of well-known pleiotropic loci for autoimmune rheumatic diseases are the HLA region, DNASEL13, TNIP1, and IRF5, among others. The identification of the causal molecular mechanisms behind the genetic associations is still a challenge. However, recent advances have been achieved through mouse models and functional studies of the loci. Here, we provide an updated overview of the genetic architecture underlying these four autoimmune rheumatic diseases, with a special focus on the HLA region.

Differential diagnosis of systemic lupus erythematosus and Sjögren's syndrome using machine learning and multi-omics data.

Systemic lupus erythematosus and primary Sjogren's syndrome are complex systemic autoimmune diseases that are often misdiagnosed. In this article, we demonstrate the potential of machine learning to perform differential diagnosis of these similar pathologies using gene expression and methylation data from 651 individuals. Furthermore, we analyzed the impact of the heterogeneity of these diseases on the performance of the predictive models, discovering that patients assigned to a specific molecular cluster are misclassified more often and affect to the overall performance of the predictive models. In addition, we found that the samples characterized by a high interferon activity are the ones predicted with more accuracy, followed by the samples with high inflammatory activity. Finally, we identified a group of biomarkers that improve the predictions compared to using the whole data and we validated them with external studies from other tissues and technological platforms.

A genome-wide association study suggests new susceptibility loci for primary antiphospholipid syndrome.

Objectives: Primary antiphospholipid syndrome (PAPS) is a rare autoimmune disease characterized by the presence of antiphospholipid antibodies and the occurrence of thrombotic events and pregnancy complications. Our study aimed to identify novel genetic susceptibility loci associated with PAPS. Methods: We performed a genome-wide association study comprising 5,485 individuals (482 affected individuals) of European ancestry. Significant and suggestive independent variants from a meta-analysis of approximately 7 million variants were evaluated for functional and biological process enrichment. The genetic risk variability for PAPS in different populations was also assessed. Hierarchical clustering, Mahalanobis distance, and Dirichlet Process Mixtures with uncertainty clustering methods were used to assess genetic similarities between PAPS and other immune-mediated diseases. Results: We revealed genetic associations with PAPS in a regulatory locus within the HLA class II region near HLA-DRA and in STAT4 with a genome-wide level of significance. 34 additional suggestive genetic susceptibility loci for PAPS were also identified. The disease risk allele in the HLA class II locus is associated with overexpression of HLA-DRB6 , HLA-DRB9 , HLA-DPB2 , HLA-DQA2 and HLA-DQB2 , and is independent of the association between PAPS and HLA-DRB1*1302 . Functional analyses highlighted immune and nervous system related pathways in PAPS-associated loci. The comparison with other immune-mediated diseases revealed a close genetic relatedness to neuromyelitis optica, systemic sclerosis, and Sjögren's syndrome, suggesting colocalized causal variations close to STAT4 , TNPO3 , and BLK . Conclusions: This study represents a comprehensive large-scale genetic analysis for PAPS and provides new insights into the genetic basis and pathophysiology of this rare disease.

Protocol for large scale whole blood immune monitoring by mass cytometry and Cyto Quality Pipeline.

Mass cytometry (MC) is a powerful large-scale immune monitoring technology. To maximize MC data quality, we present a protocol for whole blood analysis together with an R package, Cyto Quality Pipeline (CytoQP), which minimizes the experimental artifacts and batch effects to ensure data reproducibility. We describe the steps to stimulate, fix, and freeze blood samples before acquisition to make them suitable for retrospective studies. We then detail the use of barcoding and reference samples to facilitate multicenter and multi-batch experiments. For complete details on the use and execution of this protocol, please refer to Rybakowska et al. (2021a) and (2021b).

Whole blood DNA methylation analysis reveals respiratory environmental traits involved in COVID-19 severity following SARS-CoV-2 infection.

SARS-CoV-2 infection can cause an inflammatory syndrome (COVID-19) leading, in many cases, to bilateral pneumonia, severe dyspnea, and in ~5% of these, death. DNA methylation is known to play an important role in the regulation of the immune processes behind COVID-19 progression, however it has not been studied in depth. In this study, we aim to evaluate the implication of DNA methylation in COVID-19 progression by means of a genome-wide DNA methylation analysis combined with DNA genotyping. The results reveal the existence of epigenomic regulation of functional pathways associated with COVID-19 progression and mediated by genetic loci. We find an environmental trait-related signature that discriminates mild from severe cases and regulates, among other cytokines, IL-6 expression via the transcription factor CEBP. The analyses suggest that an interaction between environmental contribution, genetics, and epigenetics might be playing a role in triggering the cytokine storm described in the most severe cases.

Scoring personalized molecular portraits identify Systemic Lupus Erythematosus subtypes and predict individualized drug responses, symptomatology and disease progression.

OBJECTIVES: Systemic Lupus Erythematosus is a complex autoimmune disease that leads to significant worsening of quality of life and mortality. Flares appear unpredictably during the disease course and therapies used are often only partially effective. These challenges are mainly due to the molecular heterogeneity of the disease, and in this context, personalized medicine-based approaches offer major promise. With this work we intended to advance in that direction by developing MyPROSLE, an omic-based analytical workflow for measuring the molecular portrait of individual patients to support clinicians in their therapeutic decisions. METHODS: Immunological gene-modules were used to represent the transcriptome of the patients. A dysregulation score for each gene-module was calculated at the patient level based on averaged z-scores. Almost 6100 Lupus and 750 healthy samples were used to analyze the association among dysregulation scores, clinical manifestations, prognosis, flare and remission events and response to Tabalumab. Machine learning-based classification models were built to predict around 100 different clinical parameters based on personalized dysregulation scores. RESULTS: MyPROSLE allows to molecularly summarize patients in 206 gene-modules, clustered into nine main lupus signatures. The combination of these modules revealed highly differentiated pathological mechanisms. We found that the dysregulation of certain gene-modules is strongly associated with specific clinical manifestations, the occurrence of relapses or the presence of long-term remission and drug response. Therefore, MyPROSLE may be used to accurately predict these clinical outcomes. CONCLUSIONS: MyPROSLE (https://myprosle.genyo.es) allows molecular characterization of individual Lupus patients and it extracts key molecular information to support more precise therapeutic decisions.

Machine Learning for the Identification of a Common Signature for Anti-SSA/Ro 60 Antibody Expression Across Autoimmune Diseases.

OBJECTIVE: Anti-Ro autoantibodies are among the most frequently detected extractable nuclear antigen autoantibodies, mainly associated with primary Sjögren's syndrome (SS), systemic lupus erythematosus (SLE), and undifferentiated connective tissue disease (UCTD). This study was undertaken to determine if there is a common signature for all patients expressing anti-Ro 60 autoantibodies regardless of their disease phenotype. METHODS: Using high-throughput multiomics data collected from the cross-sectional cohort in the PRECISE Systemic Autoimmune Diseases (PRECISESADS) study Innovative Medicines Initiative (IMI) project (genetic, epigenomic, and transcriptomic data, combined with flow cytometry data, multiplexed cytokines, classic serology, and clinical data), we used machine learning to assess the integrated molecular profiling of 520 anti-Ro 60+ patients compared to 511 anti-Ro 60- patients with primary SS, patients with SLE, and patients with UCTD, and 279 healthy controls. RESULTS: The selected clinical features for RNA-Seq, DNA methylation, and genome-wide association study data allowed for a clear distinction between anti-Ro 60+ and anti-Ro 60- patients. The different features selected using machine learning from the anti-Ro 60+ patients constituted specific signatures when compared to anti-Ro 60- patients and healthy controls. Remarkably, the transcript Z score of 3 genes (ATP10A, MX1, and PARP14), presenting with overexpression associated with hypomethylation and genetic variation and independently identified using the Boruta algorithm, was clearly higher in anti-Ro 60+ patients compared to anti-Ro 60- patients regardless of disease type. Our findings demonstrated that these signatures, enriched in interferon-stimulated genes, were also found in anti-Ro 60+ patients with rheumatoid arthritis and those with systemic sclerosis and remained stable over time and were not affected by treatment. CONCLUSION: Anti-Ro 60+ patients present with a specific inflammatory signature regardless of their disease type, suggesting that a dual therapeutic approach targeting both Ro-associated RNAs and anti-Ro 60 autoantibodies should be considered.